Annihilation of vortex dipoles in an Oblate Bose-Einstein Condensate 
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We theoretically explore the annihilation of vortex dipoles, generated when an obstacle moves 
through an oblate Bose-Einstein condensate, and the possible reasons for the annihilation. We show 
that the grey soliton, which results from the vortex dipole annihilation, is lower in energy than 
the vortex dipole. We also investigate the annihilation events numerically and observe that the 
annihilation occurs only when the vortex dipole overtakes the obstacle. Furthermore, we find that 
the noise reduces the probability of annihilation events. This may explain the lack of annihilation 
events in experimental realizations. 
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PACS numbers: 03.75.Mn, 03.75.Kk 



I. INTRODUCTION 

One of the important developments in recent exper- 
iments on atomic Bose-Einstein condensates (BECs) is 
the creation of vortices and the study of their dynamics 
[U [2] . Equally important is the very recent experimen- 
tal observation of a vortex dipole, which consists of a 
vortex- antivortex pair, when an obstacle moves through 
a BEG [3] and in situ observation of vortex dipoles pro- 
duced through phase imprinting [4 . In superfluids, the 
vortices carry quantized angular momenta and are the 
topological defects, which often serve as the conclusive 
evidence of superfluidity. In a vortex dipole, vortices of 
opposite charges cancel each other's angular momentum 
and thus carry only linear momentum. This is the cause 
of several exotic phenomena like leap frogging, snake in- 
stability [5 , orbital motion [6 , trapping [7 , and others. 
The effects of vortices are widespread in classical fluid 
flow ^ and optical manipulation [9 . A good description 
of vortices in superfluids is given in Ref. [10] and review 
articles [HI [12] . More detailed discussion of vortices is 
given in Ref. [13]. 

Among the important phenomena associated with the 
Bose-Einstein condensate (BEG), the creation, dynam- 
ics, and annihilation of vortex dipoles carry lots of in- 
formation associated with the system. Several methods 
have been suggested to nucleate vortices. Recently, nu- 
cleation of the vortices have been observed experimen- 
tally by passing a Gaussian obstacle through the BEG 
with a speed greater than some critical speed [3 . The 
trajectories of these vortex dipoles are ring- structured as 
described in Refs. [MlllS]- Annihilation of vortices in 
the BEG has been speculated in a number of studies. 
However, no extensive study has been done so far. More- 
over, it has not been observed in the experiments done 
with an oblate BEG. The study of vortex dipole annihila- 
tion will shed light on the process that leads to minimum 
separation between vortex-antivortex and conditions for 
annihilation along with other phenomena arising due to 
dynamics of vortex dipoles. 

In this article, we present analytical and numerical re- 
sults related to vortex dipole annihilation for an oblate 



BEG at zero temperature. The results have been ob- 
tained using Gross-Pitaevskii (GP) equation. In the Sec- 
tion. II of the paper, we provide a brief description of the 
two-dimensional (2D) GP equation and vortex dipole so- 
lutions. Approximate analytic solutions of 2D GP equa- 
tion with a vortex dipole or a grey soliton, expressed as 
linear combinations of the eigenstates of trapping poten- 
tial, are given in Section [lllj Stability of the solutions are 
analysed from the energies of the solutions. The numeri- 
cal results, confirming the analytic results, are discussed 



in Section IV, and we then conclude. 



II. SUPERFLUID VORTEX DIPOLE 

In the mean-field approximation, the dynamics of the 
dilute BEG is very well described by the GP equation 



ihdt^iv.t) = [H + /7|^(r,t)|2]^(r,t), 



(1) 



where ^, and U are normalized wave function of 
the condensate, single-particle Hamiltonian and interac- 
tion strength respectively. The complete single-particle 
Hamiltonian H consists of a kinetic-energy operator, an 
axis-symmetric harmonic trapping potential, and a Gaus- 
sian obstacle potential, i.e.. 
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V'+--^{x'+aY+l3-'z')+Voks{x, y, t), (2) 



where a and /3 are the anisotropics along y and z axis 
respectively, and Vohs{x^ y^t) is the repulsive Gaussian 
obstacle potential. Experimentally, obstacle potential is 
produced by a laser beam which is blue-detuned with 
respect to the frequency of the atomic transition and can 
be written as 



Vohs{x, y,t) = Vo{t)ex.p 



{x — vtY + y'^ 



(3) 



where Vo{t) is the potential at the center of the Gaussian 
obstacle at time t, v is the velocity of the obstacle along 
X-axis, and wq is the width of the obstacle potential. 
In the present work, we consider the motion of obstacle 
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along X-axis only. Defining the oscillator length of the 
trapping potential aosc = \/^/ (muj)^ and considering huj 
as the unit of energy, we can then rewrite the equations 
in dimensionless form with transformations f = r/aosc, 
i = tuj, and the transformed order parameter 



(4) 



where N is the number of atoms in the condensate. 
For the sake of notational simplicity, hereafter we de- 
note the scaled quantities without tilde in the rest of the 
manuscript. In the pancake-shaped traps (a = 1 and 
/3 1)^ the order parameter 

(/)(r, t) = ^{x, y, t)C{t) exp(-z/3t/2), (5) 

where C = (/3/(27r))^/'^ exp(-/3^V4) is the order param- 
eter in the axial direction. The Eq. ([T]) can then be re- 
duced to the two dimensional form 

h Vohs[x, y,t) 



1 fd^ 

2 \dx^ 



dy'^ 
-^u\^{v,t)\^ 



d 



V^(r,t) = 0, (6) 



here u = 2aA^v^27r^/aosc, where a is the 5- wave scatter- 
ing length, is the modified interaction strength. In the 
present work, we consider condensate consisting of ^^Rb 
atoms in F = 1, mi? = 1 state with 5- wave scattering 
length equal to 99ao [16 . We have neglected a constant 
term corresponding to energy along axial direction as it 
only shifts the energies and chemical potentials by a con- 
stant without affecting the dynamics. This equation can 
be solved numerically by using Crank-Nicholson method 



A. Generation of vortex-dipole pairs 

There are several theoretical proposals to generate vor- 
tices in non-rotating traps. These include stirring of 
the condensate using blue-detuned laser or several laser 
beams [3^, H!! , adiabatic passage [18 , Raman transitions 
in bicondensate systems [19 , laser beam vortex guiding 
[20] , laser beam diffraction on a helical light grating [21] , 
and phase imprinting [22 . Among these methods, the 
easiest one to nucleate vortex dipoles is by stirring a BEC 
with a blue-detuned laser beam. When the velocity of the 
laser beam exceeds a critical velocity, vortex- antivortex 
pairs are released from the localized dip in the number 
density created due to the laser beam. These vortex 
dipoles then move through the BEC and exhibit vari- 
ous interesting dynamics [H [Ml [23] . The critical velocity 
depends on the number density, width and intensity of 
the laser beam, and the frequency of the trapping po- 
tential. This nucleation process exhibits a high degree 
of coherence and stability, allowing us to map out the 
annihilation of the dipoles. In an axis symmetric trap, 
a vortex dipole is a metastable state of superfiuid fiow 
with long lifetime. 



B. Trajectory of a vortex dipole 

The motion of a vortex dipole in a trapped BEC may 
be understood in terms of two contributions to the ve- 
locity of each vortex in the dipole. First, the precession 
due to the inhomogeneity of the condensate and secondly, 
the velocity induced by the other vortex. In our case, we 
have considered highly oblate BEC with a few well sep- 
arated vortex dipoles, hence the contribution to velocity 
field due to the other vortex dipoles can be neglected. 
Moreover, as the obstacle is moved along x-axis, it leads 
to creation of the vortex dipoles located symmetrically 
about this axis. The velocity of the each vortex in the 
dipole is then given by 

V = cj^k X r + (7) 
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where ujp is the precession frequency, k is the direction 
of circulation at the location of vortex, and 2y is the 
separation between vortex and antivortex. This equation 
can be reduced to two equations in x?/-plane as 

dx 1 



dt 
dy 
dt 



2y 



■ UJnX. 



(9) 



These two equations govern the motion of a vortex 
dipole. The solutions of these equations with different 
initial conditions are shown in Fig. [l] From the Fig. [l] 
one can observe that the vortex dipole come closer only 
along the x-axis of the BEC. They can even be closer 
than the coherence length. For this, the vortex needs to 
be created close to the obstacle. Hence, the annihilation 
will occur only along the diameter of BEC. 




X (arbitrary units) 



FIG. 1. The trajectories of the vortex dipole in oblate BEC 
calculated numerically from the equation of motion Eq. (8|9 ). 
The trajectories depend on the initial conditions. 



III. VORTEX DIPOLE ANNIHILATION 

To analyse the vortex dipole annihilation, we consider 
a model system where the vortex and antivortex are 
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static. However, we vary the distance of separation and 
examine the energy of the total system. 

The vortex and anti-vortex annihilation occurs when 
the separation is less than ^. A simplistic model of 
a vortex dipoles in the BEC of trapped dilute atomic 
gases is the superposition of harmonic oscillator eigen- 
states. The minimalist wave function which supports a 
vortex and antivortex at cordinates (— a/c, —^/hjd) and 
{—a/c^ \/h/d) is 



y) = [ia — b-\- ixc -\- dy ) 



2A 



(10) 



where a, 6, c, and / are positive variational parame- 
ters. The wave function is a superposition of the scaled 
ground state and the first and the second excited states of 
harmonic oscillator along the x and ?/-axes, respectively. 
The wave function is ideal for weakly interacting conden- 
sates, however, the qualitative descriptions remain unal- 
tered for strongly interacting condensates. The normal- 
ization condition, considering correction for the excited 
states being negligible, is + 6^ = 2/(/7r). This is a 
constraint equation. 

With a slight modification, the trajectory of the vortex 
dipole can be represented in terms of the time-dependent 
parameters. For example, the wave function when the 
vortex dipole is moving along x-axis 

ij{x, y, t) = e-'^' [-h{t) + ic{x - vt) + dy'^] e-(^'+^')//, 

(11) 

where v is the velocity with which the vortex dipole 
moves along the x-axis, and /i is the chemical potential 
of the system. 



A. Diametric vortex dipole 

Consider that the vortex and antivortex are located 
on the diameter of the condensate. Without loss of gen- 
erality, we consider the diameter as coinciding with the 
?/-axis, which is equivalent to a = in Eq. [lOj Such an 
assumption does not modify the qualitative descriptions, 
but expressions are far less complicated. The wave func- 
tion is then 



i^{x,y,t) 



-b{t) 



ICX - 



(12) 



The nontrivial phase of the wave functio n is discontinu- 
ous along X = line for — ^Jh/d ^y ^ ^Jh|d. Across the 
discontinuity, there is a phase change from — tt to tt as 
we traverse along x-axis from 0~ to 0+. To assign phase 
and elucidate the nature of the phase discontinuity, con- 
sider ^ = ^, a line parallel to x-axis, where (5 < ^Jh|d 
is a very small positive constant. Along the line y — 
[dy^ - h(i)) < and at x > 0, cx/{dy'^ - b{t)) -oo 
and is ~ 7r/2. When we decrease x, increases and it 
is ^ TT when x ^ 0+. On the other side of the ?/-axis, at 
X <C 0, cx/{dy'^ — b{t)) oo and 6 is ^ — 7r/2. However, 
on increasing x the phase tends to — tt. So, there is a 
discontinuity across the y-ax.is. This is the typical phase 



pattern associated with vortex dipoles. For the present 
case, the ground state wave function is 



iljg{x,y,t) 



-i^^+yn/f 



(13) 



the normalization condition 
J^oo I^oo li^gl'^dxdy = 1, we get the constraint equation 



and from 

-co oo 



2 

7^' 



(14) 



For general considerations, rewrite the additional term 
as 



dilj{x, y, t) = e-'^' {icxe-'"^' + dy^e'"^') 



(15) 

So that the total wave function = ?/^g + ^?/^, where 5il) 
represents an elementary excitation of the condensate. 
From the Bogoliubov theory, it satisfies the normalization 
condition 



/ [c^x^ - dPy^) e-^^^"+y"^'Uxdy = 1, 

-OO J — oo 



which after evaluation is 



(2c" - M^f) = 1. 



(16) 



(17) 



The limits, in reality, must be the physical size of the 
condensate clouds. However, due to the exponential de- 
cay and relatively large extent of the condensate clouds 
there is little difference in the results. The Eq. (17), like 



Eq. (14), is an additional constraint equation to the wave 



function. We can calculate the total energy of the sys- 
tem, without the obstacle potential, as 



/oo pO 
/ 
-OO ^oo 

+u\tp{x,y)\^ 



dxdy. 



y 



(18) 



This is the energy of condensate with a vortex dipole with 
the assumption that it is a weakly interacting system. 
Energy without the vortex may be calculated trivially 
[To]. In general, the energy added to the system due 
to the vortex dipole is not large compared to the total 
and for obvious reason, the angular momentum of the 
condensate is still zero. 



B. Grey soliton 

A slight modification to the wave function can describe 
a solitonic solution along ?/-axis. The form of the modifed 
wave function is 

V^(x, y) = [b{t) + icx + dy^] e-(^'+^')/'^, (19) 

where except for the change in the sign of 6, all the terms 
remain unaltered as in Eq. (10). It is a grey soliton as 
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the density n oc {b -\- dy'^)'^ + (a + cx)^ has a dip but is 
different from zero. 

Like in vortex dipole, we examine the nature of the 
phase around the points (0, —^/hjd) and (0, y^^bjd). The 
nontrivial phase of the wave function has no d iscontin- 
uinities. Consider again the hue y = S <^ which is 

parahel to x-axis. Along the hue {dy'^ + > at ah 
values of x and — 7r/2 ^ ^ 7t/2. In the x ^ domain, 
cx/{dy'^ + b{t)) oo and 6 ^ 7r/2. When we decrease x, 
decreases and it is when x ^ 0+. On the other side 
of the ?/-axis, at x <C the phase ~ — 7r/2. However, on 
decreasing x the phase tends to 0. So, the phase varies 
smoothly from — 7r/2 to 7r/2 along the normal to the line 
which connects (0, —y^bjd) and (0, \/h/d). 

Using the wave function in Eq. ([l9|, we can then eval- 
uate the total energy of the system E^^. We define the 
energy difference between the two states of the conden- 
sates as 



140 /im/s. Increasing the speed of obstacle generates two 
pairs in 140 /im/s ^ v < 160 /im/s and more than two 
when V ^ 160 /im/s. In other words, the number of vor- 
tex dipoles created can be controlled with the speed of 
the obstacle. Creation of vortex dipoles above a critical 
speed Vc is natural as the vortex nucleation must sat- 
isfy the Landau criterion [24]. The density and phase of 
the condensate after the nucleation of vortex dipole for 
V = 160 /im/s is shown in Fig. [sj The figure clearly 
shows nucleation dynamics of the vortex dipoles. 

Through a series of calculations, we determine Vc ~ 
90 /im/s. This is, however, less than the acoustic veloc- 
ity of the medium g, whic h depends on the local conden- 
sate density s = ^JnU jm. This also explains the reason 
for the predominant vortex dipole nucleation around the 
edge of the condensate where the density is lower and the 
acoustic speed is accordingly lower. 



AE = Evd-^gs, (20) 
which after evaluation is 

= [646^^ + Xhd^fu + 8/(8 + ^u)\ . (21) 

256 ^ 

The most general solution is when all the constants are 
positive, then l\E > and the grey soliton is lower in 
energy. This shows that when the vortex- ant ivortex col- 
lides, it is energetically favourable for them to decay into 
grey soliton. As discussed in the results section, this is 
confirmed in the numerical calculations. 

The analysis so far is for an ideal system at zero tem- 
perature, where we have neglected the thermal fluctua- 
tions and perturbations from imperfections. In addition, 
there is dissipation from three body collision losses in the 
condensates of dilute atomic gases. 

IV. NUMERICAL RESULTS 

The specific paremeters we consider for the numerical 
studies are: the species is chosen as ^^Rb with = 2 x 
10^ atoms. The trapping potential and obstacle laser 
potential parameters are same as those considered in Ref. 
[3], i.e., cj/(27r) = 8 Hz, /3 = 11.25, Vo{0) = 93.0 hw and 
wq = 10 jam. The obstacle potential Vohs is initially 
located at —12.5 aosc- To nucleate vortices, we move the 
obstacle along the x direction at a constant velocity with 
decreasing intensity until Vohs vanishes at 5.18 aosc- 



B. Vortex dipole annihilation 

It is observed that the vortex dipole annihilation is 
critically dependent on the initial conditions of the nu- 
cleation. For this reason, the annihilation events are ob- 
served only for specific range of As an example, the 
annihilation event when v is 120 /im/s is shown in Fig. [2j 
In Fig. |2j we can notice the density minima arising from 
the annihilation and propagating away from the l^bs- 

One observation, which is common to all the vortex 
dipoles getting annihilated is the nature of their trajec- 
tory. All of them traverse through Kbs, whereas the ones 
which do not get annihilated avoid Vohs- This is again 
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A. Vortex dipole nucleation 

We study the nucleation of vortices by V^bs with the 
translation speed v ranging from 80 /im/s to 200 /im/s. 
Vortices are not nucleated when the speed is 80 /im/s. 
However, a vortex-antivortex or a vortex dipole is nu- 
cleated when the speed is in the range 90 /im/s < v < 



FIG. 2. A vortex dipole is nucleated when the obstacle 
potential traverses the condensate at a speed of 120 /im/s. 
The vortex dipole, however, passes through and overtakes the 
obstacle. Later, as seen in (e), the vortex dipole annihilates 
and generates a grey soliton. The figures in the left panel 
show the density distribution and those on the right show the 
phase pattern of the condensate. 
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related to the initial conditions. The vortex dipoles are 
generally nucleated at the aft region of the Vohs where 
there is a trailling superflow. When nucleated very close 
to each other and with high velocity, the mutual force 
further increases the velocity of the vortex dipoles. At 
the same time, it decreases the distance separating vortex 
and antivortex. So, the kinetic energy is high enough to 
surpass lobs- Later, at some point vortex and antivortex 
separation is less than ^, and they annihilate. 
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FIG. 3. A vortex dipole is nucleated as the obstacle potential 
traverses the BEG with a speed of 160 /xm/s. The figures 
in the left panel shows the density with time, where time 
progresses from top to bottom. Figures on the right panel 
show the phase pattern of the condensate. 



A reliable and qualitative way to describe occurrence 
of annihilation could be achieved by observing the den- 
sity at the cores of vortex and antivortex which form the 



dipole. For the vortex, the matter density at the core 
when V is 120 /im/s is shown in Fig [4j In the plot, at 
Time ~ 900 (arbitrary units), the core density starts in- 
creasing. This is because the core starts to fill in with 
the atoms from around the vortex after the annihilation. 
As discussed in Section. 



annihilation can occur 
only when it is energetically favorable. In other words, 
the state with the grey soliton must have lower energy 
than the vortex dipole. This is clearly seen in the change 
of total energy of the system, which is shown in Fig. [Sj 
After the annihilation, the energy of the system decreases 
and continues to do so at a steady rate. Although not 
shown in the plot, before the annihilation the energy is 
on an average constant. The energy change, can thus be 
considered as a very good signature of the annihilation 
of vortex dipoles. 




80.8902' 

100 200 300 400 500 600 700 800 900 1000 

Time (arbitrary units) 

FIG. 5. Variation in the energy of the system with time (in 
scaled units). There is a decrease in energy after the vortex 
dipole annihilation, and the plot is based on the results of cal- 
culations in which a vortex dipole is imprinted at (—2.0 ttosc, 
±0.1 aosc)- It is then allowed to annihilate and evolve in time. 
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FIG. 4. Density variation at the core of the vortex with time 
(in scaled unit). After the vortex dipole annihilation, density 
increases till it reaches the bulk value. The values correspond 
to 120 /im/s obstacle speed. 



It is to be mentioned that for the parameters consid- 
ered in the present work, the speed of sound is 2190 /im/s, 
which is similar to the speed observed by Neely et. al. 
[3]. The coherence length of the system is ~ 0.25 /im and 
the minimum distance between the vortex dipole pair is 
~ 0.3 /im. The position of annihilation determined from 
the equation of motion (Eq. [8][9| matches with those ob- 
tained from the simulation. 



C. Effect of noise and dissipation 

In the numerical studies, the annihilation events are 
not rare. But, this is in contradiction with the exper- 
imental results of Neely and collaborators [3 , they ob- 
served no signatures of annihilation events. One possible 
reason is that our numerical calculations are too ideal, 
and one immediate remedy is to include fluctuations. For 
this we introduce white noise during the real time evo- 
lution. One immediate outcome is, the symmetry in the 
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trajectory of the vortex and antivortex is lost. The su- 
perflow around the vortex is no longer a mirror reflection 
of the antivortex, which was nearly the case without the 
white noise. The deviations are shown for an example 
case in the Fig. |6] The change in path leads to the 
suppression of annihilation of vortex dipoles. 





(b) 

• 




f(c) 





X (a 
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to the Hamiltonian 1-L. Based on the previous work [25] 
for ^^Rb, the inelastic two-body loss rate coefficient 
K2 = 4.5 X 10~^^ cm^ s~^, and inelastic three-body loss 
rate coefficient Ks = 3.8 x 10~^^ cm^ s~^. With trap 
loss, the annihilation events continue to occur. However, 
during the destructive time of flight observations in the 
experiments, the decreased atom numbers may lower 
the contrast and reduce the possibility of observing an 
annihilation event. 



V. CONCLUSIONS 



max 



2k 



FIG. 6. The trajectory of a vortex dipole in the presence of 
white noise. There is a lack of symmetry in the trajectory 
of the vortex and antivortex. This reduces the possibility 
of an annihilation event significantly. The figures in the left 
(right) panel show the density (phase) of the condensate and 
time increases from top to bottom figures of each panel. The 
speed of the obstacle is 180 /xm/s, and the white noise is at 
the level of 0.01%. 



The other important effect is the loss of atoms from the 
trap. We have examined the effect of loss terms, which 
arise from inelastic collisions in the condensate. There 
are two types of inelastic collisions that lead to the loss 
of atoms from the trap: two body inelastic collision loss 
and the three body loss. To model the effect of loss of 
atoms from the trap, we add the loss terms 



■i^3|*(r,t)r 



(22) 



When an obstacle moves through a condensates above 
a critical speed, it nucleates vortex dipoles and the num- 
ber of dipoles seeded depends on the obstacle velocity. 
Depending on the initial condition of nucleation, vortex 
and antivortex annihilation events occur under ideal con- 
ditions: at zero temperature, no loss, and without noise. 
However, the noise destroys the superflow reflection sym- 
metry around the vortex and antivortex. This reduces 
the possibility of annihilation events and may explain the 
lack of annihilation events in experimental observation in 
Ref. [3]. 
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